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We propose a simple geometrical construction of topological invariants of 3-strand Brownian 
braids viewed as world lines of 3 particles performing independent Brownian motions in the complex 
plane z. Our construction is based on the properties of conformal maps of doubly-punctured plane 
z to the universal covering surface. The special attention is paid to the case of indistinguishable 
, particles. Our method of conformal maps allows us to investigate the statistical properties of the 

' topological complexity of a bunch of 3-strand Brownian braids and to compute the expectation 

. value of the irreducible braid length in the non-Abelian case. 
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I. INTRODUCTION 



Our paper presents some results concerning geometry and statistics of random three-strand braids. Besides the 
relevance of this work at a purely algebraic level (see, for example, for a review of principal topological problems), 
, its usefulness in the understanding statistical physics of entangled lines, independent on their physical nature, is 
I also undeniable. Statistics of ensembles of uncrossible linear objects with topological constraints has very broad 
' application area ranging from problems of self-diffusion of directed polymer chains in flows and nematic-like textures 
to dynamical and topological aspects of vortex glasses in high-Tc superconductors 0,01 • In order to have representative 
and physically clear image for the system of fluctuating lines with non-Abelian (i.e. noncommutative) topology we 
formulate the model in terms of entangled Brownian trajectories: such representation serves also as a geometrically 
clear image of Wilson loops in (2-|-l)D non-Abelian field-theoretic path integral formalism. In particular our paper 
O , focuses on the geometrical, algebraic and statistical properties of the simplest non-commutative braid group B^. 
^ . 

Consider the ensemble of N particles randomly moving in the complex plane z — x + iy. Let us label the coordinates 
of these particles by Zj (t) — Xj (t) + iyj (t) where 1 < j < N and t is the current " time" (the initial configuration of 
the particles corresponds to t = 0). It is clear that in (2+l)D "space-time" (z, t) the diffusive motion of all particles is 
ci ' described by the statistics of N "world lines" or "directed polymers" and the time t plays the role of the length of the 
"world line". In what follows we shall assume the periodic boundary conditions, i.e. we shall suppose that at some 
time t = T the configuration of the particles in the plane z is just the permutation of the initial configuration, i.e. 
{zi(T), Z7v(T')} = 7'{zi(0), ZAr(O)}, V being a permutation of N elements. Thus, under the natural condition 
Zj{t) 7^ Zk{t) for any 1 < {fc,j} < N, the "world lines" shall be entangled. 

The topology of the braid can play a crucial role in macroscopic physical properties. Let us mention only two 
examples. The elastic properties of polymer networks strongly depend on the initial degree of entanglement between 
subchains forming the sample and, hence, the elastic properties of the rubbers can be controlled by different initial 
conditions of samples preparation 0, Q- In Cu02-based high-Tc superconductors in fields less than the critical 
magnetic field Hc2 there exists a region where the Abrikosov flux lattice is molten, but the sample of the superconductor 
demonstrates the absence of the conductivity. This effect is explained by highly entangled state of flux lines due to 
their topological constraints 0, Q . 

We restrict our study to = 3, i.e. to the case of 3-strand braids produced by bunches of "word lines" of 3 particles 
simultaneously moving in the complex plane z (see fig^. The simplest non-commutative group is defined on the 
projection of this bunch onto the plane (t, x) as it is shown in fig^ The group B^ is constituted by the set of 
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generators {cti, 0-2, CTj^ ^, (Tj"^} satisfying commutation relations explicited below. Each topological configuration can 
be exactly encoded by an element of the braid group B^, that is a "word" written in terms of "letters" -the braid group 
generators, associated to elementary moves-positive/negative crossings in the projection (see figQ). A configuration 
of the braid corresponds at the time t = T to an element ujt of -B3. 

We are interested in an explicit construction of topological invariants of entanglements of such world lines i.e., in the 
3~strand random "braid" lot- We investigate and evaluate the complexity of such randomly generated braid, defined 
as the minimal number of generators L(lot) necessary to write lot- This quantity L{ujt) is called the irreducible 
length (in the metric of words) and exactly coincides with the minimal number of crossings necessary to represent the 
entanglement. Therefore L{lot) can be chosen as an indicator of the braid complexity. In particular, if L{ujt) — the 
braid is trivial (i.e. the strands are unentangled). 




FIG. 1: Braiding of three directed lines. Elementary moves ai in the projection generate the group B3. 

In what follows we shall repeatedly use the so-called Burau matrix representation of the braid group. In the 
particular case of the group B3 the Burau representation is given by 2 x 2 matrices: 

where i is a free parameter. It is known that for B3 this representation is faithful When t = —1 the group B3 
coincides with the modular group PSL{2,Z). 

The paper is structured as follows. The first part introduces the basic notions of homotopy groups and homotopic 
length, necessary to self-contained description of topological invariants. The method proposed in this part, demon- 
strates on the simplest examples how the monodromy representations of groups allows one to build the universal 
coverings, giving rise to topological invariants, associated at a physical level to fluxes of non-Abelian extension of 
"solenoidal magnetic fields". Using the same ideas, the second part deals with the case of B3 for which we explicitly 
derive a topological invariant, directly linked to the irreducible braid length, and construct its non-Abelian flat con- 
nection. For the model of ideal Brownian braids, the stochastic evolution of this topological invariant is given and its 
asymptotic distribution is derived. 
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II. TOPOLOGICAL INVARIANTS AND MONODROMY REPRESENTATIONS 



A. Basic concepts and definitions 



Let us consider the double-punctured complex plane C** = C — {zi, Z2}, of variable z = x + iy and suppose the 
coordinates of the punctures Mi and M2 to be Zi = (0, 0) and Z2 = (c, 0) correspondingly. We set c = 1 by appropriate 
rescaling of the plane z. 

Take now two closed elementary paths on z such that the first one (71) encloses only the point Mi and the second 
one (72) surrounds only the point M2. With the usual composition law of paths, 71 and 72, generate a fundamental 
group F2, which is the first homotopy group, tti, of the double-punctured complex plane z: 



The trivial path (the unit element of the group F2) is the composition of an arbitrary loop with its inverse: e = 
7i7j~^ = 7j~^7i, i = {1,2}. The loops 7^ and 7^ are equivalent if 7i can be continuously deformed into 7i (the 
equivalent loops represent the same element of ^2)- We denote the class of equivalent paths by 7^. 

Any element of the group F2 being finitely generated, corresponds to a closed path on z and can be represented 
by a "word" consisting of a sequence of letters-generators of the group F2: {71, 72, 7r^, 72^^}- Each word can be 
reduced to a minimal (or "irreducible") representation. For example, the word W = 7172^^71 71 7r^72'^727r^72~^ = 
7172"^ [71 [71 7r^] [72'^72] 7r^] 72"^ can be reduced to = 7i72'^72'^- Consider now a finitely generated group 
G = (51,52) (to be more precise, the group admitting only one- or two-generators presentation). We are interested 
in the monodromy representation of F2 into G defined by the following group homomorphism 



Note that this straightforwardly can be generalized to the case of multi-punctured surfaces or other topological spaces. 
We do not discuss in details the existence of such homomorphism, assuming that it results mainly from the fact that 
F2 is a free group. 

Our main goal in the present paper consists in defining an explicit geometrical construction of a non-Abelian 
generalization of a Gauss topological invariant (i.e. a " linking number" ) for different groups G using the monodromy 
representations of G. In particular, we construct a complex flat connection A{z) on C**, i.e. a "Bohm-Aharonov-like 
vector potential" A{z), whose holonomy gives rise to a representation of G. A special attention is paid to the case 
when G = Bs- Moreover, we use the developed approach to estimate the averaged complexity of a 3-strand braid 
represented by independent motion of 3 Brownian particles in (2+1) dimensions. 



In order to explain the basic notions, we begin with the simplest possible case — the Abelian group Z and construct 
the corresponding topological invariant by means of conformal transforms. 



7ri(e*) = J^2 




(2) 



B. Topological invariants from conformal maps 



1. Abelian case: commutative group G = Z 



The central point of the approach deals with the reconstruction of a linear differential equation on the manifold C** 
by its monodromy group G. In other words, defining the action of G on the space of solutions of some 2'^'^-order liner 
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differential equation with two branching points, we are attempting to recover the form of this differential equation 
and its solutions. In the case of G = Z it is natural to consider a usual 1-dimensional additive action 

n[w(z)] = w{z) + 2inn 

for w defined on C** and n G Z. It turns out that the simplest linear differential equation with two branching points, 
satisfying those conditions is: 

(Z- Zl)(z- Z2)-^ - 2Z + Zl +Z2 = (3) 

whose solution reads 

= ln(z — Zi) + ln(z — z^) (4) 

One can check that the free group r2 acts on the space of solutions as follows 

w(z) '^TiijJ-^ w{z) + 2i7r (5) 

which means that the homomorphism ^ is trivial in for this example: ^'(7i) = \ (i = 1,2). The function vj{z) 
conformally maps the doubly punctured plane to the universal covering space w — u^iv free of any branching points. 
In the complex plane w we have 

= In |(z - zi)(z - Z2)| 
: arg(z - zi) + arg(z - Z2) 

The function u'(z|zi, Z2) is a topological invariant because of equality In particular, the function v{z\z\^ z-i) defines 
the total number of turns in the plane z around the branching points z\ and Z2 and hence it is nothing else as the 
Gauss linking number. 

Thus, knowing the conformal transform w(z) of the multiply punctured plane C** to the universal (i.e. uniformizing) 
covering surface, we can easily extract a topological invariant Inv(7) of a closed path 7 (starting and ending at some 
arbitrary point zq ^ {0, 1} in the plane z) from the difference Win(zo) — Wfin(-Zo)- The path 7 connects the images of 
the point zq on different Riemann sheets — the "copies" of the fundamental domain of G. Recall that by definition 
the fundamental domain of a group G is a minimal connected domain tessellating the whole covering space under the 
action of G. Representing the topological invariant Inv(^) (7) as a full derivative along the contour 7, we get: 

Inv(2)(7) =^ itifin - Win = ^ '^^dz' ^ / A{z)dz (6) 

7 

The physical interpretation of the derivative A{z) = ^^^^^ is very straightforward. The conformal transform w{z) 
plays the role of a complex potential of a field A{z), which defines a flat connection of a multiple-punctured plane 

For the commutative group G — Z we obtain A{z) by taking the derivative of Q): 

Aiz) = + (7) 

Z — Zl Z — Z2 

This expression can be easily identified with the standard ( Abelian) Bohm-Aharonov vector potential of two solenoidal 
magnetic fields orthogonal to the plane z and crossing it in the points zi and Z2. 

In the next section the same construction of the flat connection associated to a specific group G will be generalized 
to the non-Abelian case. 
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2. Non-Abelian case: non- commutative groups F2, Hq 



We consider a special class of hyperbolic and hyperbolic-like groups: the free (F2) and the Hecke {Hq) groups, 
as well as the braid group B^. By hyperbolic-like groups we mean a class broader than hyperbolic groups in the 
classification of M.Gromov (According to the Gromov's definition, the group does not belong to the class of 
hyperbolic groups). The important feature for us would be just the exponential growth of the group. From this point 
of view the group fits our scheme. 

1. The free group F2 (in general, Fn (n = 2,3,4, ...)) by definition is the free product of two (in general of n) copies 
of cyclic groups of second order, Z2. The matrix representation of the generators of the free groups F2{/i,/2} and 
F3{h,S} are: 




(8) 



2. The Hecke group Hq is the free product of two cyclic groups Z2, and Zq of orders 2 and q respectively. The Hecke 
group is defined by the relations 

iSTqr=b^q^l 

where the generators Tq and S have the following matrix representation 

/ 1 2 cos ^ \ / 1 \ 

The parameter q takes discrete values q = 3,4, 5,6, . . .. The Hecke group Hq "interpolates" between the modular 
group PSL{2,Z) (for q = 3) and the free group F3 with 3 generators {q 00). 

We have stressed in the previous section that the topological invariant can be constructed on the basis of the 
conformal map w(z) of multiple-punctured plane to the universal covering space of a group. We now extend the 
described method to more interesting cases than considered at length of the section Hi fj II 

We derive the conformal mapping of the half-plane Im z > onto the fundamental domain of the triangular group 
G — a curvilinear triangle lying in the upper half-plane Imw > 0. The action of G on this fundamental domain 
generates the whole covering space. Each copy of the fundamental domain represents a Riemann sheet corresponding 
to the fibre bundle above z and the whole covering space w is the unification of all such Riemann sheets — see figlH 

The coordinates of initial and final points of a trajectory on the universal covering w determine: 
• The coordinates of the corresponding points on z; 



The element of G corresponding to the homotopy class of the path on z. 



Following the construction described in the previous section let us define the action of G{gi,g2} in the complex 
plane w = u + iv. The group G admits the faithful 2-dimensional representations and acts in the covering space w by 
fractional-linear transforms. Consider two basic contours 71 and 72, associated to the action oi G on z {z ^ {zi, Z2, 00}. 
The contours 71 and 72 enclose the branching points located at zi and Z2 correspondingly (z 7^ {zi,Z2,oo}). The 
function w{z) {z 7^ {zi, Z2,oo}) obeys the following transformations: 



w 



71 



. aiw{z) + bi r 72 1 ~ , . a2w{z)+b2 

wi(z) — — ; w z^z W2{z) — -— — (11) 

^ ^ Ciw{z)+di' I J ^ C2w{z)+d2 ^ ' 
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FIG. 2: Conformal transforms of the upper half-plane Imz > to the fundamental domain of the group F2 (up) and of the 
group P5L(2,Z) (down). 

The matrices gi and (72 

9i=("'^']; 32=f"'^0 (12) 
\ci di J \C2 d2 J 

are the matrices of basic substitutions of the group G{gi,g2}, i.e. gi and g2 are the generators of G. 

It is well known that the function w{z) can be defined as the quotient of two fundamental solutions ui{z) and U2{z) 
of a second order differential equation with branching points {zi — (0,0), Z2 = (0, 1),2;3 — 00}. As it follows from 
the analytic theory of differential equations, the solutions ui{z) and U2{z) undergo the linear substitutions when the 
variable z moves along the contours 71 and 72: 

yu2{z) J yu2{z) J yu2(z) j yu2{z) j 

where 171 and 32 are the generators of the monodromy group of this equation. The problem of restoring the form of 
a differential equation by the monodromy matrices gi and g2 of the group G of the differential equation, is known as 
the Riemann-Hilbert problem. Yet we restrict ourselves with the groups F2 and Hq, the group B3 shall be considered 
separately. 



The free group F2. For the free group F2 the solution of the Riemann-Hilbert problem gives rise to the following 
second-order differential equation |l4|: 

, d'^u(z) , , duiz) 1 , , , , 

Indeed, a possible basis of solutions of this equation is as follows: 

?/i(.2) = F(l/2, 1/2,1, z) 

(15) 

U2(z) =iF(l/2,l/2,l,l-z) 
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92 = 1 





Using the well-known properties of hypergeometric functions, one can restore the monodromy matrices defined in 
for this basis: 

(16) 

which coincides with the generating set of the group i^2- The function w{z) performing the conformal map of the 
upper half-plane Imz > onto the fundamental domain (the curvilinear triangle ABC) of the universal covering w 
satisfies ea.()ll|l and can be written as: 

w{z) = —-^ 17) 



The function w{z) is well known in the literature (see, for example, 'l4l) and its inverse z{w) is the elliptic modular 
function 



z{w) ^ = (18) 



Now we can give an explicit expression of the flat connection A(z) for the doubly punctured plane corresponding to 
the monodromy of the free group F2- Taking the derivatives and using the properties of the hypergeometric functions, 
we get: 

_ dw{zl E{z)K{l -z) + (E{1 - z) - K{\ - z))K{z) 
dz ~ 2z{z-l)K^{l~z) 

where 

7r/2 ir/2 

K{z) = f : E(z) = [ deVi-zsin^e (20) 

^ 

are correspondingly the complete elliptic integrals K{z) and E{z) (see, for example, jl5j). 

In the FigOlwe have plotted the absolute values |^(2)| of the Abelian (Eq.Q) and the non-Abelian fEa.l|19()') flat 
connections. 
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Hence even in the vicinity of the branching point z = the function A{z) defined in (|19|l does not coincide asymp- 
totically with the flat connection for the Abelian case iQ. 

The Hecke group Hq. Wc now pass to the case of Hq, keeping in mind that for q ~ 3 we recover the group 
PSL{2, Z), directly linked to B3 {B^ is a central extension of PSL{2, ZV see The solution of the corresponding 

Riemann-Hilbert problem leads to the following differential equation: [ig 



z{z~l) 



d u{z, q) 



dz^ 

whose possible basis of solution is: 



2 qj \ q 



du{z, q) 1 



dz A\2 q 



- --- u{z,q)^Q (22) 



(23) 



where 



r(i - i)r(a±2)r(3g+2N 
^ ^ r(i + i)r(23^)r(2£z2) ^ ^ 



Taking into account that the series defining the hypergeometric functions F(z) converges for \z\ < 1 and corresponds 
to the so-called logarithmic case (see |3) one obtains from (|23|l the following monodromy matrices after proper 
analytic continuation: 

One can check that after normalization by the determinant, the matrices ()25|l become the generators of Hq (corre- 
spondingly of orders q and 2). Let us point out that for q — 3, one can identify the inverse function z(w) with the 
Klein's "absolute modular invariant" 

Z[W) _ J[W) - 5408(0^ ^)^8(o,^)0|(O,u;) 

where 62, 03, are the Jacobi elliptic ^-functions. As in the previous section, the corresponding complex flat connec- 
tion Afig {z) can be obtained as the full derivative 

. . X dwH, jz) 

from the conformal map 

ui{z) 

U2{Z) 

with ^1(2;) and U2{z) defined in (|23|l . 



III. INVARIANT OF 3-STRAND BRAIDS 



Let us return to the model of random braids discussed in the very beginning of the paper. The system of three 
independently moving particles is described by 3 complex variables (i.e. 6 degrees of freedom) and hence our system 
can not be directly viewed as a monodromy problem. In order to describe the whole system by one complex variable C 
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and to be able to use the tools elaborated in the previous section, with the minimal loss of information, we introduce 
the anharmonic quotient: 

^23 ft) _ Z2it) - Z3{t) 

One can easily check that the variable ^ contains the complete topological information of a mutual configuration of 
three entangled lines, except the global phase, or global twist of the braid. It means that we do not take into account 
the center of B3 considering the factor group PSL{2,Z) = B^fE. We can always take into account the global twist 
afterwards by passing to the rotating coordinate system with the origin located in the center of mass of the system 
of three particles. From the statistical point of view the possibility of neglecting the center of B3 (which is equivalent 
to considering the factor group PSL{2,'Z) only), has been discussed in the papers 0,^3- iB,ct, in these paper it 
has been shown that the escape rates for random walks on B3 and PSL(2,'L) are the same in the limit of infinitely 
long trajectories. 

Thus, the parametrization of the three-particle system {zi(t), Z2(i), z^{t)} by the function C(t) living in the doubly 
punctured complex plane enables us to preserve all topological characteristics of the braid of mutually entangled world 
lines of these particles. However the precise derivation of the expression for the flat connection, as it is shown below, 
explicitly depends on the fact whether the particles are identical or not. 

Looking at the elementary moves associated with the generators ui, 02 of -B3, we obtain the transformations for 
the variable For example, let us pick up some path in the homotopy class of cti (see fig^l. By definition of cti, we 
have at some time t: 



ai[zi ^ Z2, Z3 ^ Z3] : < 



^ z-i{t) = 2:3 

In the same way we can get the transformation of C along a path corresponding to the homotopy class 02 = 02 \z\ ^ 
zi,Z2 <-> 2:3]. Finally we arrive at the following set of transformations of Q: 

^ 1 

which precisely define the representation of the symmetric group S'3. 



Indistinguishable particles. Our main goal consists in constructing the flat connection for the system of three 
identical particles moving in the plane. The condition of the indistinguishability of particles requires to factor the 
action of the group B3 by S'3, and hence to consider C as living in the factor space C**/S'3. The case of distinguishable 
particles shall be discussed at the end of this section. 

All the trajectories ({t) (parameterized by t) are obviously closed in the space C**/S3. Keeping in mind our strategy 
of the previous section, we want to define a monodromy representation of B3 (more precisely of PSL{2, Z)) acting on 
C** /S3. To do that we construct a conformal map Ws(C) of the doubly punctured plane factorized over the action of 
the symmetric group, C**/S'3 onto the fundamental domain of the modular group, PSL(2,7j): 

WsiO : eV5'3 ^ PSL{2,Z) (29) 



To our knowledge, the explicit method of solving the uniformization problem H29|) has not been yet considered in 
the context of our problem. So, our own way to tackle it deals with the following observation. 
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Consider the so-called pure braid group P3 defined as follows: 

P3 = {alai,{a,a2f) (30) 

It is obvious that the group P3 is the subgroup of B3. In particular, the group P3 can be identified with the factor 
group F2 X 1j through the obvious correspondence 

= 91, = 51 (31) 

where gi and (72 are the generators of the free group F2. The relations H14(l follow from the Burau representation of 
braid group generators 1^ and the obvious geometric construction shown in the fig^ 




12 3 12 



FIG. 4: The squares of the braid group generators at act as the free group generators gt: gi = (i — 1, 2). 

In the case of the pure braid group P3 there is no difference between distinguishable and indistinguishable particles 
moving in the plane z because the generators gi and 52 correspond to " full turns" of one world line of the particle 
with respect to another one (and not to "half turns" as in case of the braid group B3). That is, we deal with the 
monodromy representation of ^2- Let us remind that we factor out the global turns (i.e. the center of the group) in 
the space C** . The uniformization problem in this case is solved by the conformal map Ws (C) of the doubly punctured 
plane to the covering space 

WsiO : F2 

Thus, we return to the topic of the previous section, where an explicit form of such conformal map Ws(C) = w{(^) is 
given in (|17|l . 

Consider now the whole group P3/Z, take the same function w{() and check the effect of the basic substitutions of 
the symmetric group (|28ll onto the transformations of w{C). Using the well known properties of the modular fmictions 
(see, for example, |2(i |21i, (2^ ) , we get: 



( C ^ 



(32) 



technically it is much easier to check (instead of (|32|l ') the transformations for the inverse function Ci'^s) ^ ^^(o.ma) . 



C-IJ w,^l 

The transformations of Ws precisely coincide with the basic transformations of the modular group PSL(2, Z). In fact, 

Ws + 1) ^ i 




ij C-1 

Hence the function Ws{C) explicitly solves the desired uniformization problem H29II . Thus, we can use this function 
:{C) to construct the non-Abclian flat connection As{C) for the trajectories on C** /S3 parametrized by ((t): 



MO = ^ (34) 
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where 



U2iO F(l/2,l/2,l,l-z) 
(see eq.lO). 

The expression H34(> is the source for the exphcit construction of the topological invariant via equation of three 
entangled world lines zi(t), Z2{t) and Z3(i) shown in fig^ 



Distinguishable particles. For distinguishable particles the action of the symmetric group 5*3 should be neglected 
and hence C is living in the doubly punctured plane, C**. The problem of uniformization in this case is solved by the 
conformal map ^^(C) mapping the doubly punctured plane, C**, onto the fundamental domain of the modular group, 
PSL{2,Z): 

WdiC) : C** — > PSL{2,Z) (35) 

(compare this expression to (|29|l 'l. The corresponding function ^^(C) coincides with W2{C) — ^1[q found in the 
privious section (use the equations (|23|) - (|24|l for q = 3). 



IV. STOCHASTIC BEHAVIOR OF THE INVARIANT FOR BROWNIAN 3-STRAND BRAIDS 



We are interested in this section in the stochastic behavior of the invariant |w(?7)| when each of the three particles 
perform a 2D Brownian motion (BM) represented by its complex coordinate Zi. This model naturally appears as 
the generalization of the Edwards' problem of determining the distribution of the winding angle of two independent 
2D BM The intrinsic noncommutative structure of the 3-particles model, contrary to the commutative 2- 
particles problem, induces a completely different behavior of the invariant, which characterizes quantitatively the 
tangle complexity. The method employed here is inspired by a work of Gruet [2^ . 



A. The Ideal Model 



We begin with an ideal case of point-like particles diffusing in the infinite plane with the diffusion constants D set 
to 1. Following previous sections, one has to study the motion of the reduced variable ry in the space C**/S'3. Our 
strategy relies on the following theorem on conformal transforms of complex BM often referred to as a generalized 
Paul Levy's theorem j^]. Denote by Bt the generic BM in the complex plane z. For any conformal map /, the image 
Yt = f{Bt) of a BM Bt is a time changed BM Bu{t) such that 

du = ^dYdY = \f'{Bt)\^dt (36) 

In terms of the diffusion equation (which corresponds to the Langevin equation H3t)|l . the time change can be interpreted 
as a space dependent diffusion coefficient |25l |: 

dii 

dtP{Yt) = ^mYt) (37) 

A slightly modified version of the generalized Levy theorem shows that being the ratio of two complex BM 

(see (|26|) '). is also a time changed BM, denoted ({t) = Bu{t)- The time change then reads, following the rules of Ito 
calculus: 

du — -dCdC — ^ i^J dt (38) 
4 \zz2{t)V 
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Now we should pass to the covering space w and describe the evolution of the invariant w{() = w{Bu)- Using 
once more the generalized Levy's theorem, we obtain that there exists a hyperbolic BM B such that w(C) is again a 
time-changed B: 

wiCit)) = (39) 

with 

... K(M_,„ 

Im-(to(B,.)) 

The choice of the hyperbolic metric is not artificial and anticipates the geometrical properties of the transform w(C)- 
Note the importance of obtained result: we know that the topological invariant w{Q, giving the braid complexity, 
performs a hyperbolic diffusion in a new time, or equivalently a hyperbolic diffusion in a space with the metric- 
dependent diffusion coefficient. We now have to study the full time change T{t). Combining (|38|) and (|40|l and using 
the well known inverse function = z{w) — X{w) (see (IIHI) ). we obtain after separating the variables the following 
implicit relation between t and r, which defines the functional Ht{T): 

H,{r) = f = r ^ltlLW(B.,)rfr' (41) 

This expression is exact, but does not allow a straightforward interpretation. Let us first notice that the i-dependence 
can be easily shown to be, with the probability one: 

lim ^ Ci > (42) 

t^oo hit ^ ' 

The sensitive point is the r-dependence. The integral kernel / 

is invariant under the transformations of the group F2, and Htir) is therefore an additive functional of the BM on the 
quotient surface EI/F2. This motion is ergodic. Therefore if Ht{T) is integrable over the fundamental domain of F2, 
one can clear up the asymptotic r-dependence (see in particular a similar approach proposed by Gruet in 26] ). In 
our case Ht (r) is not integrable over the domain of F2 ; the r-dependence for this ideal model requires more attention 
and have to be treated with more care. We expect to consider that question in details separately, while below we 
express a conjecture based on the fact that the non-integrability does not lead to a growth faster than r. We indeed 
believe that the non-integrability gives rise to sub-leading terms in r and the following limit holds: 

hm = C2 (44) 

Comparing (|42|) and (|44|1 we arrive at the relation 

li. . §i (45) 

t^oo t Lj2 

Using (|45|l we formulate a central result connecting the hyperbolic BM B^ describing the connection of the expec- 
tation of our topological invariant d{w) being the distance in the covering space w, with the word length L{Bt) of 
the corresponding random braid. Recall that the distribution of B^. is well known, and in particular, 

lim^ = l 

r—^oo T 

almost surely. Moreover, if one notices that the word lengths L for PSL{2. Z) and for F2 are quasi- isometric, we have: 



(HBr)) 

lim — — — C3 > in probability 

T^oo rlnr 



(46) 
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what is a direct consequence of a theorem by Gruet proved in |2g . Combining this result with we end up with 
the following leading asymptotics at large time for the expectation of the word length of Brownian 3-strand braids: 

(Lt) - ln(i) X In(lni) for < ^ c» (47) 

The average topological complexity increases with time, however the growth is very slow. It is limited (as in the 
case of two particles where this logarithmic scaling also holds) by the angular part of the 2D BM in an infinite space 
|2{t1 I27I 1 : when two particles are far apart (and the corresponding typical distance grows with time) , their relative 
angle varies very slowly. 



B. Model with a compact domain 



We here discuss the physical reasons that justify a regularization of the integral 141|l over the fundamental domain 
of F2. The divergence of this integral occurs when ij approaches the points {0, 1, 00} meaning that either two particles 
collide, or one particle goes to infinity. It is then quite natural to modify slightly the model in order to avoid these 
pathological cases, which in fact are not realistic if one describes the physical objects of finite thickness such as 
polymers. Specifically, we assume that three Brownian particles evolve in a bounded domain of characteristic size R 
and that they experience hard-core repulsion at distance r. This can be re-expressed in terms of constraints due to 
hard walls added to the domain of rj, hereafter denoted by D: 

e{r,R) < \t]\ < M{r,R) 

(48) 

|??(r,i?)-l| >e 

The precise shape of the boundary of D is not important. Note that 'w{D) is compact in H. Restricting our model 
to this domain, we can now claim the integrability of the functional HP{t) over the domain w{D). Within this 
approximation, we obtain the following r-dependence from the ergodic property offered above: 

lim ^^i^ = C4 > (49) 

r — >oc T 

In this model the i-dependence is readily changed as follows: 



Comparing H49|) and (|5U|I we arrive at 



ttD ( \ 

lim = C5 > (50) 

t^oo t 



lim -Y^T^ (51) 



The expectation of d{w) can now be described asymptotically: d(w) is a hyperbolic BM in the variable t oc t (see 
(ISlfl ). In particular, we have with the probability one: 

lim ^ > (52) 

t^oo t 

For a compact domain one can straightforwardly check that the topological invariant, d, and the word length, L, are 
quasi-isometric. The asymptotic behavior of the expectation value of L is therefore similar to the one of {d(w{ri))): 

(if) - i for i ^ 00 (53) 

For the bounded domain, the topological complexity grows much faster, as the infinite space effect discussed above is 
absent. The scaling (if) c>c t should be compared to the (commutative) scaling (icomm) < x. \/ i for two particles. It 
is noteworthy to stress that H53|) is fully consistent with a discrete model considered in [l8|, [lOj . 
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V. 



CONCLUSION 



To conclude, we would like to comment the new physical content of our results. So far the only topological invariant 
studied in the context of entanglements of fluctuating linear objects was the so-called winding number ili, id] j 
which is known to be incomplete for more than two linear objects. We here propose a model, which describes exactly 
the underlying non-Abelian topology of the problem. 

We propose a simple geometrical construction of topological invariants of 3-strand Brownian braids viewed as world 
lines of 3 particles performing independent Brownian motions in the complex plane z. Our construction is based on the 
properties of conformal maps of doubly-punctured plane z to the universal covering surface. We pay special attention 
to the case of indistinguishable particles. Our approach is mainly "self made" and its geometrical transparency we 
consider as the basic advantage. The standard machinery of constructing the nontrivial braid group representations 
from the Conformal Field Theory is outlined in the Appendix where we mainly review the strategy realized in [2^ . 
The Appendix is added to our paper exclusively to establish some links between the approachs based on the geometry 
of conformal maps and on CFT. 

Our method of conformal maps allow us to investigate the statistical properties of the topological complexity of a 
bunch of 3-strand Brownian braids and to compute the expectation value of the irreducible braid length. 



In this section we try to describe the derivation of the braid group representation from the monodromics of some 
CFT. We will mainly describe the system under investigation in physical terms, rather than in more rigorous but less 
transparent algebraic topological setting. 

So, consider n (n = 3 shall be mainly treated) indistinguishable particles, living in the complex plane C. The "static" 
(or configurational) phase space is the following set X„ = YnjSn = CxF„_i/S'„ where F„ = {z = (zi, z„) S C"; i ^ 
i Zi 7^ Zj \ and S'„ is the group of permutations of n elements. The above defined decomposition means that one can 
factor out the coordinate of the center of mass of the system. Thus, we study the topology of Ar„, or more precisely 
of its first homotopy group 7ri(A"„), describing the classes of equivalence of closed curves on A"„ where i3„ = 7ri(X„). 

In physical terms we can rephrase the said above as follows. A quantum mechanical description of this system 
requires to define the wave functions, or in other words an n fold tensor-valued functions ^' : X„ y®" over this 
configuration space. These functions should represent the internal quantum numbers such as a spin for each particles. 
The most general case consists in taking as the so-called G-modules, where G is a Lie group. The obtained 
multi-valued (in fact, tensor-valued) structure is called a fiber bundle. The main difficulty in such construction is as 
follows: the functions ^'(z) and 5'(z' 7^ z) belong to different spaces, say and /i^', and therefore can not be directly 
compared. To overcome this difficulty one has to define a map , '■ ^ h^i that "transports" the wave function 
from z to z' along the path 7. 

Now let us define the holonomy operator, ,, and a one-form lo = dTzz, called "flat connection" over X„. The 
important point is that the flatness of the connection 



is a necessary and sufficient condition for the holonomy group at z to give rise to a monodromy representation of 
the fundamental group 7ri(Ar„). Recall that the holonomy group is the group of T^^^ of all closed paths with the 
natural path composition as an internal law. We are therefore lead by this statement to a study of flat connections 
over Xn- 

What is the physical meaning to be extracted from such consideration? The basic idea, put forward by the discovery 
of the non-Abelian Bohm-Aharnov effect |3l(, is that a closed trajectory in the phase space can affect internal quantum 




APPENDIX A: BRAID GROUP REPRESENTATION AND CFT 



cLu + Lu A uj — 
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numbers of the system giving rise to topological interaction which depend only on the homotopy class of this closed 
path. The flat connections usually considered in the literature are the following matrix valued 1-forms, or the so-called 
Knizhnik-Zamolodchikov (KZ) connections: 



.KZ _ 

\<<a<h<n 



W^^= 2^ 17„fcdl0gZab (Al) 



where ilab € End(y®") is the Casimir invariant (not explicitly written here), and Zah = Za — Zb- The reader is referred 
to li^l for a more detailed description of this aspect. This connection was introduced in 2D conformal field theory in 
the context of the Wess-Zumino-Novikov-Witten model and is related to the study of chiral current algebras [sj [s^ . 
In this model, a primary field is covariant under two kinds of transformations: local gauge transformations generated 
by the current J and conformal reparameterization generated by the stress-energy tensor T. These T and J are linked 
by the Sugawara formula describing consistency between two covariances, and leads hence to the KZ equation for an 
n-points correlator ^': 

hd^ = uj-i' (A2) 

h being a complex parameter. This formulation suggests a geometric interpretation: 'if is covariantly constant for 
this connection. The equation ljA2|l being a first order linear differential equation, can be treated in the frameworks 
of the analytic theory of differential equations. Any solution of (|A2p can be represented by a linear combination of 
its 71 — 1 fundamental solutions. The formal (and in practice not very useful) description of the holonomy group: 

T^^^ = Pexp / uj^^ (A3) 

then reduces to a matrix representation. It is shown in that equation l)A2|l can be reduced in the case n = 3 to a 
system of ordinary differential equations for the S'[/(A^)-invariant amplitude F I 7] = — ) defined by 

V ^13/ 

*(zi, Z2, ^3) = ^"3* (?7(1 - v))'^F{v) (A4) 



h^ + J^l^-^]Fi7,) = (A5) 



This parameterization yields 

dr] 1 — T] rj 

where Ki2^ K23 G End(F®"). Denoting then by /q, Ii the basis of SU{N) invariant tensors in V^^ such that 

K23I0 - 0, X12/1 = (A6) 

one can write 

F{v) = (1 - v)fiv)Io + Vf\v)h (A7) 
and then reduce (|A2p to the following system of ordinary differential equations [2^: 



arj 
arj 



(A8) 



which admits the form of ordinary 2nd order Riemann differential equation with branching points at 77 = 0, 1 
d'P ^((. 4\ 2 \ dp 



ri{v-l)^^+{U-T)-^-^ + lr^^J^+(l-^\[l-'i]P=Q (, = 0,1) (A9) 
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The basis of fundamental solutions of (|A9|) can be written in terms of standard hypergeometric functions 

^"('^^ = ""m- ' ' vr^'+^-'FU - h^\j + h^\j + 2h-\ 1 ~ V) 

_i ' _i (AlO) 

mv) = ''^%h-\2h'-!'P '^'''^''^^' ' i-j+h-\i-j+ 2h-\^) 



Using (|A10p and the integral representations of the hypergeometric functions, the authors of |23| have directly compute 
the action of B3 generators cti, (T2 in the basis of solutions. Choosing a path in the homotopy class of ai {i = 1, 2): 

Z' ~\~ Z' 1 

ai : z^.it) = — — + 2'^«.'+i'^~"*' = const (All) 

corresponding to an elementary move Ui, they end up with the following braid relations: 

a. : fl - mkfi (A12) 

where the monodromy matrices are: 

with q = . This is a 2-dimensional representation of B3. 

Following the general method developed at the length of the Sections im and UTTI we can conjecture that 



is a topological invariant of the group B3. In particular, we expect that — - — gives the corresponding explicit 

drj 

expression of the non-Abclian flat connection. 



[1] C. Belisle, Ann. Prob. 4, 1377 (1989) 

[2] A. Comtet, J. Desbois, C. Monthus, J. Stat. Phys. 73, 433 (1993) 

[3] D.R. Nelson, Phys. Rev. Lett. 60, 1973 (1988); D.R. Nelson, H. Seung, Phys. Rev. B 39, 9153 (1989) 
[4] D.R. Nelson, A. Stern, Proceedings of the XIV Sitges Conference, in Complex Behavior of Glassy Systems, June 10-14 
(1996) 

[5] M. Doi and S.F. Edwards Theory of Polymer Dynamics (Academic Press: NY, 1986) 
[6] B. Drossel, M. Kardar, Phys. Rev. E 53, 5861 (1996) 

[7] L. Treloar, The Physics of Rubber Elasticity, 3rd ed., (Clarendon Press: Oxford, 1975); M. Gottlieb, R.J. Gaylord, Polymer, 
24, 1644 (1983); A.R. Khokhlov, F.F. Ternovskii, Sov. Phys. JETP, 63, 728 (1986) 

[8] S. Obukhov, M. Rubinstein, Phys. Rev. Lett. 65, 1279 (1990) 

[9] Epstein, Word Processing in Groups (Barlett and sons: Boston, 1992) 
[10] S. Nechaev, A. Grosberg, A. Vershik, J. Phys. A, 29, 2411 (1996) 
[11] P.Dehornoy, Braids and self-distributivity (Birkhauser, 2000) 

[12] J. Birman, Braids, Links and Mapping Class Groups (Ann. Math. Studies 82, Princeton Univ. Press, 1976) 
[13] M. Gromov, Hyperbolic groups (Springer: New York, 1987) 

[14] Higher Transcendental Functions eds.: H. Bateman, A. Erdelyi (McGraw Hill, 1955) 

[15] Handbook of mathematical functions, eds.: M. Abramowitz and I. Stegun (Dover Publications, Inc.: New York, 1970) 

[16] A.B. Venkov, Preprint LOMI P/2/86 

[17] N. Ikeda, H. Matsumoto, J. Funct. Ann., 163, 63 (1999) 

[18] R. Voituriez, Nucl. Phys. B 621, 675 (2002) 



17 



[19] S. Nechaev, R. Voituricz, J. Phys. A, 36, 43 (2003) 

[20] V.V. Golubev Lectures on Analytic Theory of Differential Equations (Moscow: GITTL, 1950) 
[21] E. Hille, Analytic Function Theory (Ginn & Co.: Boston, 1962) 
[22] J. Harnad, J. McKay, Proc. Roy. Soc. Lond. 456 261 (2000) 
[23] D.R. Hofstadter, Phys. Rev. B, 14, 2239 (1976) 
[24] A. Comtet, Ann. Phys., 173, 185 (1987) 

[25] K. Ito, H.P. McKcan, Diffusion Processes and their Sample Paths (Springer: Berlin, 1974) 
[26] J.C. Gmct, Prob. Thcor. Rcl. Fields 111, 489 (1998) 
[27] S.F. Edwards, Proc. Phys. Soc. London 91, 513 (1967) 

[28] D. Revuz, M. Yor, Continuous Martingales and Brownian Motion (Springer: Berlin, 2001) 
[29] I. Todorov, L. Hadjiivanov, Phys. Atom. Nucl. 64, 2059 (2001) 
[30] A. Tsuchiya, Y. Kanie, Lett. Math. Phys., 13 303 (1987) 

[31] E. Verlinde, A note on braid statistics and the nonabelian Aharonov-Bohm effect, in Modem Quantum Field Theory, 

(World Scientic: Singapore, 1991) 
[32] E. Witten, Coomm. Math. Phys., 92, 455 (1094) 
[33] D.G. Knizhnik, A.B. Zamoodchikov, Nucl. Phys., B247, 83 (1984) 



